3.1990 \(\int \sqrt{a+\frac{b}{x^3}} x^2 \, dx\)

Optimal. Leaf size=47 \[ \frac{1}{3} x^3 \sqrt{a+\frac{b}{x^3}}+\frac{b \tanh ^{-1}\left (\frac{\sqrt{a+\frac{b}{x^3}}}{\sqrt{a}}\right )}{3 \sqrt{a}} \]

[Out]

(Sqrt[a + b/x^3]*x^3)/3 + (b*ArcTanh[Sqrt[a + b/x^3]/Sqrt[a]])/(3*Sqrt[a])

_______________________________________________________________________________________

Rubi [A]  time = 0.0828673, antiderivative size = 47, normalized size of antiderivative = 1., number of steps used = 4, number of rules used = 4, integrand size = 15, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.267 \[ \frac{1}{3} x^3 \sqrt{a+\frac{b}{x^3}}+\frac{b \tanh ^{-1}\left (\frac{\sqrt{a+\frac{b}{x^3}}}{\sqrt{a}}\right )}{3 \sqrt{a}} \]

Antiderivative was successfully verified.

[In]  Int[Sqrt[a + b/x^3]*x^2,x]

[Out]

(Sqrt[a + b/x^3]*x^3)/3 + (b*ArcTanh[Sqrt[a + b/x^3]/Sqrt[a]])/(3*Sqrt[a])

_______________________________________________________________________________________

Rubi in Sympy [A]  time = 6.97374, size = 39, normalized size = 0.83 \[ \frac{x^{3} \sqrt{a + \frac{b}{x^{3}}}}{3} + \frac{b \operatorname{atanh}{\left (\frac{\sqrt{a + \frac{b}{x^{3}}}}{\sqrt{a}} \right )}}{3 \sqrt{a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  rubi_integrate(x**2*(a+b/x**3)**(1/2),x)

[Out]

x**3*sqrt(a + b/x**3)/3 + b*atanh(sqrt(a + b/x**3)/sqrt(a))/(3*sqrt(a))

_______________________________________________________________________________________

Mathematica [A]  time = 0.0722618, size = 67, normalized size = 1.43 \[ \frac{1}{3} x^{3/2} \sqrt{a+\frac{b}{x^3}} \left (\frac{b \tanh ^{-1}\left (\frac{\sqrt{a} x^{3/2}}{\sqrt{a x^3+b}}\right )}{\sqrt{a} \sqrt{a x^3+b}}+x^{3/2}\right ) \]

Antiderivative was successfully verified.

[In]  Integrate[Sqrt[a + b/x^3]*x^2,x]

[Out]

(Sqrt[a + b/x^3]*x^(3/2)*(x^(3/2) + (b*ArcTanh[(Sqrt[a]*x^(3/2))/Sqrt[b + a*x^3]
])/(Sqrt[a]*Sqrt[b + a*x^3])))/3

_______________________________________________________________________________________

Maple [C]  time = 0.037, size = 3340, normalized size = 71.1 \[ \text{output too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  int(x^2*(a+b/x^3)^(1/2),x)

[Out]

-1/3*((a*x^3+b)/x^3)^(1/2)*x^2/a^2*(6*I*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+
(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/
2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/
3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticF((-(I*3^(1/2)-3)*x*a/(I*
3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+
1)/(I*3^(1/2)-3))^(1/2))*3^(1/2)*x^2*a^2*b-6*I*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)
/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(
I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2
*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(1/2)-3)
*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-3),((I*
3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*3^(1/2)*x^2*a^2*b-1
2*I*(-a^2*b)^(1/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2
)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^
(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a
*x+(-a^2*b)^(1/3)))^(1/2)*EllipticF((-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^
2*b)^(1/3)))^(1/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/
2))*3^(1/2)*x*a*b+12*I*(-a^2*b)^(1/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-
a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)
+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3)
)/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(1/2)-3)*x*a/(I*3
^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-3),((I*3^(1/2)+3
)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*3^(1/2)*x*a*b+6*I*(-a^2*b)^(
2/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*
(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*
((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1
/3)))^(1/2)*EllipticF((-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(
1/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*3^(1/2)*b-
6*I*(-a^2*b)^(2/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2
)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^
(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a
*x+(-a^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a
^2*b)^(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-3),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*
3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*3^(1/2)*b-6*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(
-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*
3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b
)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticF((-(I*3^(1/2)-3)*x*
a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(
1/2)+1)/(I*3^(1/2)-3))^(1/2))*x^2*a^2*b+6*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*
x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(
1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(
1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(1/2)-3)*x*a/
(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-3),((I*3^(1/
2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*x^2*a^2*b+12*(-a^2*b)^(1
/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(
-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*(
(I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/
3)))^(1/2)*EllipticF((-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1
/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*x*a*b-12*(-
a^2*b)^(1/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*
3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3))
)^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a
^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^
(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-3),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2
)+1)/(I*3^(1/2)-3))^(1/2))*x*a*b-6*(-a^2*b)^(2/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(1/2)
-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(1/3)
)/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-
a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticF((-(I*3^(1/2)-
3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(
I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*b+6*(-a^2*b)^(2/3)*(-(I*3^(1/2)-3)*x*a/(I*3^(
1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b)^(
1/3))/(I*3^(1/2)+1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*((I*3^(1/2)*(-a^2*b)^(1/3)-2*a*
x-(-a^2*b)^(1/3))/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2)*EllipticPi((-(I*3^(
1/2)-3)*x*a/(I*3^(1/2)-1)/(-a*x+(-a^2*b)^(1/3)))^(1/2),(I*3^(1/2)-1)/(I*3^(1/2)-
3),((I*3^(1/2)+3)*(I*3^(1/2)-1)/(I*3^(1/2)+1)/(I*3^(1/2)-3))^(1/2))*b-I*(a*x^4+b
*x)^(1/2)*(1/a^2*x*(-a*x+(-a^2*b)^(1/3))*(I*3^(1/2)*(-a^2*b)^(1/3)+2*a*x+(-a^2*b
)^(1/3))*(I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3)))^(1/2)*3^(1/2)*x*a^2+3*
x*(a*x^4+b*x)^(1/2)*a^2*(1/a^2*x*(-a*x+(-a^2*b)^(1/3))*(I*3^(1/2)*(-a^2*b)^(1/3)
+2*a*x+(-a^2*b)^(1/3))*(I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3)))^(1/2))/(
x*(a*x^3+b))^(1/2)/(I*3^(1/2)-3)/(1/a^2*x*(-a*x+(-a^2*b)^(1/3))*(I*3^(1/2)*(-a^2
*b)^(1/3)+2*a*x+(-a^2*b)^(1/3))*(I*3^(1/2)*(-a^2*b)^(1/3)-2*a*x-(-a^2*b)^(1/3)))
^(1/2)

_______________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \[ \text{Exception raised: ValueError} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(a + b/x^3)*x^2,x, algorithm="maxima")

[Out]

Exception raised: ValueError

_______________________________________________________________________________________

Fricas [A]  time = 0.367058, size = 1, normalized size = 0.02 \[ \left [\frac{4 \, a x^{3} \sqrt{\frac{a x^{3} + b}{x^{3}}} + \sqrt{a} b \log \left (-{\left (8 \, a^{2} x^{6} + 8 \, a b x^{3} + b^{2}\right )} \sqrt{a} - 4 \,{\left (2 \, a^{2} x^{6} + a b x^{3}\right )} \sqrt{\frac{a x^{3} + b}{x^{3}}}\right )}{12 \, a}, \frac{2 \, a x^{3} \sqrt{\frac{a x^{3} + b}{x^{3}}} - \sqrt{-a} b \arctan \left (\frac{2 \, \sqrt{-a} x^{3} \sqrt{\frac{a x^{3} + b}{x^{3}}}}{2 \, a x^{3} + b}\right )}{6 \, a}\right ] \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(a + b/x^3)*x^2,x, algorithm="fricas")

[Out]

[1/12*(4*a*x^3*sqrt((a*x^3 + b)/x^3) + sqrt(a)*b*log(-(8*a^2*x^6 + 8*a*b*x^3 + b
^2)*sqrt(a) - 4*(2*a^2*x^6 + a*b*x^3)*sqrt((a*x^3 + b)/x^3)))/a, 1/6*(2*a*x^3*sq
rt((a*x^3 + b)/x^3) - sqrt(-a)*b*arctan(2*sqrt(-a)*x^3*sqrt((a*x^3 + b)/x^3)/(2*
a*x^3 + b)))/a]

_______________________________________________________________________________________

Sympy [A]  time = 7.51932, size = 48, normalized size = 1.02 \[ \frac{\sqrt{b} x^{\frac{3}{2}} \sqrt{\frac{a x^{3}}{b} + 1}}{3} + \frac{b \operatorname{asinh}{\left (\frac{\sqrt{a} x^{\frac{3}{2}}}{\sqrt{b}} \right )}}{3 \sqrt{a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x**2*(a+b/x**3)**(1/2),x)

[Out]

sqrt(b)*x**(3/2)*sqrt(a*x**3/b + 1)/3 + b*asinh(sqrt(a)*x**(3/2)/sqrt(b))/(3*sqr
t(a))

_______________________________________________________________________________________

GIAC/XCAS [A]  time = 0.252727, size = 53, normalized size = 1.13 \[ \frac{1}{3} \, \sqrt{a x^{4} + b x} x - \frac{b \arctan \left (\frac{\sqrt{a + \frac{b}{x^{3}}}}{\sqrt{-a}}\right )}{3 \, \sqrt{-a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(a + b/x^3)*x^2,x, algorithm="giac")

[Out]

1/3*sqrt(a*x^4 + b*x)*x - 1/3*b*arctan(sqrt(a + b/x^3)/sqrt(-a))/sqrt(-a)